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Abstract 

Axillary bud outgrowth determines shoot architecture and is under the control of endogenous hormones and a fine- 
tuned gene-expression network, which probably includes small RNAs (sRNAs). Although it is well known that sRNAs 
act broadly in plant development, our understanding about their roles in vegetative bud outgrowth remains limited. 
Moreover, the expression profiles of microRNAs (miRNAs) and their targets within axillary buds are largely unknown. 
Here, we employed sRNA next-generation sequencing as well as computational and gene-expression analysis to iden- 
tify and quantify sRNAs and their targets in vegetative axillary buds of the biofuel crop sugarcane [Saccharum spp.). 
Computational analysis allowed the identification of 26 conserved miRNA families and two putative novel miRNAs, as 
well as a number of trans-acting small interfering RNAs. sRNAs associated with transposable elements and protein- 
encoding genes were similarly represented in both inactive and developing bud libraries. Conversely, sequencing 
and quantitative reverse transcription-PCR results revealed that specific miRNAs were differentially expressed in 
developing buds, and some correlated negatively with the expression of their targets at specific stages of axillary 
bud development. For instance, the expression patterns of miR159 and its target GAMYB suggested that they may 
play roles in regulating abscisic acid-signalling pathways during sugarcane bud outgrowth. Our work reveals, for the 
first time, differences in the composition and expression profiles of diverse sRNAs and targets between inactive and 
developing vegetative buds that, together with the endogenous balance of specific hormones, may be important in 
regulating axillary bud outgrowth. 
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Introduction 

Plant architecture is denned as a three-dimensional organiza- 
tion of the plant body including shoot development pattern, 
shape, and position of leaves and flower organs (Reinhardt 



and Kuhlemeier, 2002). The aerial plant architecture is basi- 
cally determined by shoot branching and tillering, which, in 
turn, are coordinated by the formation of the shoot apical 
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meristem during the embryonic phase and subsequently by 
the formation of axillary meristems (Schmitz and Theres, 
2005; Leyser, 2009). Axillary meristems eventually develop 
into axillary buds that can arrest their growth and become 
dormant/inactive in several plant species (Ongaro et al, 
2008). After the release of shoot apical meristem dominance, 
dormant/inactive axillary buds may resume their growth 
and produce entire new axes (Shimizu-Sato and Mori, 2001; 
Leyser, 2009). Thus, axillary bud outgrowth is crucial to con- 
trol shoot architecture and biomass production. High bio- 
mass production is a desirable characteristic for biofuel crops 
such as sugarcane. 

Sugarcane (Saccharum spp) is characterized by a highly 
complex polyploid genome with publicly available expressed 
sequence tag (EST) and Genomic Survey Sequences 
Databases that assist genetic studies on this tropical crop. 
Commercial cultivars are usually propagated using stem cut- 
tings that contain three to four axillary buds. After planting, 
buds grow out to form primary shoots and develop roots to 
support the newly developing plant. Therefore, efficient bud 
outgrowth is also important in establishing and propagat- 
ing new sugarcane plantlets in the field. Notable progress 
has been done to elucidate the molecular and physiological 
bases of vegetative bud outgrowth in model plants such as 
Arabidopsis, rice, and pea (Leyser, 2009; Doust, 2007), but in 
sugarcane this developmental process is poorly characterized. 

Bud outgrowth is a complex phenomenon that relies on 
modifications in cell growth and proliferation, accompanied 
by drastic changes in the production and perception of hor- 
mones, including auxins, cytokinins, gibberellins, abscisic 
acid (ABA), and strigolactones (Leyser, 2009; Wang and Li, 
2008). Reinitiation of bud growth is determined mainly by 
the plant's genetic programme, which probably includes small 
RNAs (sRNAs). Interestingly, several members of gene fami- 
lies associated with shoot development are targets for regu- 
lation by sRNAs (Chuck et al, 2007; Koyama et al, 2007; 
Schwarz et al, 2008; Jiao et al, 2010; Wang et al, 2010). 

sRNAs (19-25 nt) are classified into two major categories: 
small interfering RNAs (siRNAs) and microRNAs (miR- 
NAs). Both are recognized as important regulators of gene 
expression and epigenetic regulation in diverse plant species 
(Voinnet, 2009; Khraiwesh et al, 2010). miRNA loci are 
transcribed by RNA polymerase II into primary transcripts, 
which are processed by nuclear RNase III -like enzyme Dicer- 
like (DCL) proteins, such as DCL1. Mature miRNAs are 
incorporated into Argonaute complexes that target degrada- 
tion or translational repression of mRNAs. As a result, miR- 
NAs play important roles in plant development, including 
cell patterning and organ development, as well as hormone 
signalling (Voinnet, 2009). 

siRNAs are a complex pool of sRNAs in plants, which 
includes frans-acting siRNAs (tasiRNAs) and repeat-associ- 
ated siRNAs (rasiRNAs) (Vaucheret, 2006; Vazquez, 2006). 
Most siRNAs target the same locus from which they were 
derived, except for phased 21 nt tasiRNAs, which are pro- 
duced by TAS loci and target mRNAs from different loci, sim- 
ilar to miRNAs (Allen et al, 2005). Phased trans- acting 21 nt 
RNAs can also be generated by other loci in plant genomes, 



probably triggered by miRNA-directed cleavage of the tran- 
scripts (Zhai et al., 2011). rasiRNAs are generated from and 
target repetitive sequences including transposable elements 
and satellite and microsatellite repeats, as well as centromere 
heterochromatic repeats (Chen et al., 2006; Domingues et al., 
2012). Recently, it has been shown that a particular class of 
transposons, miniature inverted-repeat transposable elements 
(MITEs), can generate miRNA-like rasiRNAs that have roles 
during plant stress responses and hormone signalling (Yan 
et al, 2011). 

sRNAs are often identified by computational analysis and/ 
or experimental approaches such as cloning and sequenc- 
ing (Wang et al, 2009; Wei et al, 2009; Calvino et al, 201 1). 
However, sRNA data are scarce for sugarcane, partly because 
its genome sequence is presently unavailable. As a result, few 
miRNAs have been identified in this important biofuel crop 
thus far. Recently, our research group computationally identi- 
fied 19 distinct miRNA precursors (Zanca et al, 2010), whilst 
other groups have cloned and sequenced a number of small 
RNAs associated with the drought response in sugarcane 
(Ferreira et al, 2012; Thiebaut et al, 2012). Additionally, 
at present, it is unknown whether important regulators like 
miRNAs play roles during vegetative axillary bud outgrowth. 
To enrich our knowledge of sRNAs associated with axillary 
bud development, we generated sRNA libraries from sugar- 
cane dormant/inactive and outgrown buds. We found a select 
group of sRNAs expressed in both inactive and developing 
buds, including a number of rasiRNAs and 21 nt tasiR- 
NAs. We identified numerous known microRNAs, as well 
as two novel miRNA candidates. Several miRNA targets 
were computationally predicted, and the expression of five 
was monitored during bud outgrowth by quantitative reverse 
transcription-PCR (qRT-PCR). As with their mRNA targets, 
we examined the expression profiles of selected microRNAs 
using stem-loop qRT-PCR. Further in situ hybridization 
analysis on selected miRNAs revealed distinct spatial locali- 
zation patterns within the axillary buds. Additionally, we 
quantified ABA and catabolite levels in sugarcane inactive 
and developing buds. Our work supplies novel insights into 
the dynamic developmental processes underlying vegetative 
axillary bud development, which are orchestrated at least in 
part by sRNA-mediated gene regulation. 



Materials and methods 

Plant material 

Six-month-old sugarcane plants of the hybrid SP80-3280 were used 
to obtain vegetative axillary buds. For all experiments, we excised 
axillary buds from the fourth to the eighth node of sugarcane culms, 
which showed similar outgrowth rates (data not shown). We also 
collected axillary bud tissues from mature plants of Saccharum 
officinarum (accession Muntok, Java) and Saccharum spontaneum 
(accession SES205A) grown in the field. Developing axillary buds at 
2 and 5 days after planting (DAP) were obtained from stem cuttings 
planted into an appropriate substratum and grown under green- 
house conditions. Vegetative apexes were collected from 3-week- 
old in vitro plants of the hybrid SP80-3280. Vegetative tissues of 
Sorghum bicolor (BTx623), 1-month-old plantlets of Oryza sativa 
(ssp.japonica cv. Nipponbare) and Brachypodium disthacyon (Bd21), 
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2-week-old Solanum lycopersicum (cv. Micro-Tom) seedlings, and 
2-week-old seedlings of Arabidopsis thaliana (Columbia) were also 
harvested. Arabidopsis seedlings were grown in a growth chamber 
with a 16h/8h light/dark photoperiod at 200 umol mr 2 s -1 . 

Public sequence data sets 

Sugarcane EST sequences were retrieved from the Sugarcane Gene 
Index (release 3.0, July 2011), and sugarcane bacterial artificial 
chromosome (BAC) sequences were downloaded from GenBank. 
All recorded miRNAs from plants were obtained from miRBase 
(version 17; http://www.mirbase.org/). We also employed BLAST 
searches to identify sRNAs sequences similar to tRNAs [GtRNAdb: 
http://gtrnadb.ucsc.edu/ (Chan and Lowe, 2009) and rRNAs (Rfam 
database, http://rfam.sanger.ac.uk/ Gardner et al., 201 1)]. Complete 
genomic sequences of sorghum, rice, and Arabidopsis were down- 
loaded from Phytozome (http://www.phytozome.net/), GenBank, 
and TAIR (http://www.arabidopsis.org/), respectively. Their tran- 
scriptomes were obtained from The Gene Index Project (http://com- 
pbio.dfci.harvard.edu/tgi/). 

sRNA library construction and bioinformatic analysis 
Total RNA from 0 DAP inactive and 2 DAP developing axil- 
lary buds was resolved on a 15% TBE/urea polyacrylamide gel to 
obtain the sRNA fraction of 16-30 nt. Briefly, 5'-adenylated single- 
stranded adapter was first ligated to the 3' end of the RNA using 
T4 RNA ligase without ATP, followed by a second single-stranded 
adapter ligation at the 5' end of the RNA using T4 RNA ligase in 
the presence of ATP. The resulting products were fractioned on a 
10% TBE/urea polyacrylamide gel and used for cDNA synthe- 
sis and PCR amplification. Illumina sequencing was performed 
as described previously (Blevins et ah, 2011). Raw sequences were 
retrieved in a FASTQ formatted file and the adapter sequences were 
removed using Perl Scripts. After trimming of adapter sequences, 
the inserts were sorted into separate files according to their lengths. 
Subsequently, we grouped all similar sequences and counted the 
read contributions from each library for each non-redundant sRNA 
sequence. We used the program MAQ (http://maq.sourceforge.net/) 
to map 19-25 nt RNA reads against several databases utilized in 
this study as described above. MAQ is a program that rapidly aligns 
short reads to reference sequences. The sRNA library data from 
axillary buds were deposited in the GEO database (http://www.ncbi. 
nlm.nih.gov/geo/) under accession number GSE38894. 

To identify EST sequences similar to MITEs, we firstly blasted 
the Sugarcane Gene Index sequences against the Gramineae Repeat 
Database (http://plantrepeats.plantbiology.msu.edu/index.html) using 
the following parameters: 80% of EST coverage and Expect (E) value 
<e-5. The recovering sequences were blasted against each other and 
those showing >90% nucleotide identity were regarded as one MITE- 
containing transcript. sRNAs from both libraries were blasted against 
the MITE-like EST sequences using the following parameters: no mis- 
matches allowed, word size of seven, and E value of 1000. 

tasiRNAs derived from sugarcane transcripts were identified using 
the computational approach described by Chen et al. (2007). This 
approach is non-biased and applicable to predict transcripts yield- 
ing phased sRNAs, especially when analysing data from organisms 
with incomplete genomic information, such as sugarcane. Sugarcane 
21 nt sRNAs and EST sequences were retrieved from the respec- 
tive databases as described above. The coordinates of the sRNAs 
were then analysed using the TAS prediction algorithm based on 
Perl. ESTs with P <0.001 were considered as potential tasiRNA- 
generating transcripts. 

Identification of novel and known miRNAs expressed in 
sugarcane axillary buds 

The sRNAs sequences were mapped using MAQ against miRbase 
version 17 and currently available sugarcane miRNAs sequences. To 



minimize the noise, we did not consider, for further analyses, low- 
abundance miRNAs (with a total number of reads <15). We also 
mapped the sRNAs against non-coding sequences of the Sugarcane 
Gene Index for identification of novel miRNAs in sugarcane. 
Precursor sequences of approximately 620 nt were extracted (300 nt 
upstream and 300 nt downstream from the BLAST hits) and used 
for hairpin structure predictions using the MFOLD3.2 algorithm. 
The number of structures, free energy, miRNA-like helicity, num- 
ber of arms per structure, size of helices within arms, and size/sym- 
metry of internal loops within arms were analysed by our in-house 
MIRcheck-based script (Zanca et ah, 2010), followed by manual 
inspection. To distinguish new miRNAs from rasiRNAs, we blasted 
precursor and mature sequences of sRNAs against the Gramineae 
Repeat Database (Ouyang and Buell, 2004). Finally, to identify dif- 
ferentially expressed miRNAs between 0 and 2 DAP sRNA deep- 
sequencing data, we firstly normalized the data by the number of 
the total mapped reads (transcripts per million) and then applied 
statistical significance using Fisher's exact test with Bonferroni cor- 
rection and adjusted P values <0.01 as a criteria. For target predic- 
tion, we used approaches described previously (Zanca et al, 2010; 
Dai and Zhao, 2011). 



RNA extraction and stem-loop pulsed RT-PCR 
Total RNA was extracted using Trizol reagent (Invitrogen, USA) 
according to the manufacturer's instructions and treated with 
DNAse I (Invitrogen) to remove any residue of genomic DNA. 
DNAse-treated RNA (1.5 ug) was reverse transcribed to generate 
the first-strand cDNA according to the method of Varkonyi-Gasic 
et al. (2007). Oligo(dT) primer was also added to the reaction for the 
detection of target mRNAs and internal controls. cDNA dilutions 
were used for PCRs containing 1.0 ul of cDNA, 1.5mM magne- 
sium sulfate, 0.25 mM of each dNTP, 10 pmol of each primer, and 
1 U of Tag DNA polymerase (Promega, USA). The reactions were 
done with the following cycling conditions: 94°C for 2min, and an 
appropriate number of cycles of 94 °C for 20 s, 60 °C for 30 s, and 
72 °C for 45 s. All reactions were repeated twice with two biological 
samples. Primer sequences are described in Supplementary Table SI 
at JXB online. 



Stem-loop pulsed qRT-PCR 

First-strand cDNA was transcribed as described above. A SYBR 
Green PCR was performed using a 7500/7500 fast Sequence 
Detection System (Applied Biosystems, USA). Briefly, 5 ul of 1:80 
(v/v) cDNA dilutions was added to 12.5 ul of Platinum SYBR 
Green qPCR superMix-UDG (Invitrogen), 3 pmol of each primer 
and ddH 2 0 to a final volume of 25 ul. The reactions were ampli- 
fied for 2min at 50 °C and 2min at 95 °C, followed by 40 cycles 
of 95 °C for 15 s and 60 °C for 30 s. The sugarcane 25S rRNA 
gene (TC148086) was used as an internal control (Iskandar et al., 
2004). PCR products for each primer set were subjected to melt- 
curve analysis, confirming the presence of only one peak on ther- 
mal dissociation generated by the thermal-denaturing protocol. The 
PCR efficiency for each primer set was evaluated through standard 
curve constructed from the amplification values of serial dilutions 
of cDNA. Three replicates were analysed for each biological sam- 
ple along with template-free reactions as negative controls. The 
threshold cycle (C T ) was determined automatically by the instru- 
ment, and the fold change of each gene was calculated using the 
equation 2~ AACT (Livak and Schmittgen, 2001). Primer sequences are 
described in Supplementary Table SI. 



sRNA in situ hybridization 

We used 0 and 5 DAP vegetative axillary buds harvested from 
6-month-old sugarcane hybrid SP80-3280 plants for sRNA in situ 
hybridizations as described by Javelle and Timmermans (2012). 
Locked nucleic acid (LNA) probes with sequences complementary 
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to selected miRNAs and negative controls (Scramble-miR) were 
synthesized by Exiqon (USA) and labelled with digoxigenin using 
a DIG Oligonucleotide 3'-end Labeling kit (Roche Applied Science, 
USA). Ten picomoles of each probe was used for each slide, and 
hybridization and washing steps were performed at 55 °C. 



Comparative sequence analysis 

Comparative analysis of possible orthologues of SsDCLl 
(TCI 43 136) was done by constructing a phylogenetic tree contain- 
ing highly similar plant sequences. Phylogenetic relationships of 
aligned protein sequences were constructed using the neighbour- 
joining method in MEGA4 software (Tamura et al, 2007). The 
accession numbers of the Arabidopsis and rice DCL sequences have 
been published elsewhere (Nogueira et al, 2009). 

ABA and catabolite quantification 

Vegetative axillary buds of each developmental stage (0, 2, and 5 
DAP) were harvested from 6-month-old sugarcane hybrid SP80- 
3280 plants. The procedures for tissue sample extraction, purifi- 
cation, and further quantification of ABA and catabolites were 
performed as described previously using an ultra-performance liq- 
uid chromatography tandem mass spectrometry method (Chiwocha 
et al., 2003). 



Results 

Overview ofsRNAs expressed in vegetative 
axillary buds 

Vegetative axillary bud outgrowth determines shoot architec- 
ture and is under the control of endogenous hormones and 
a fine-tuned gene -expression network. To unravel additional 
components of regulatory networks operating within axillary 
buds, we sought to identify sRNA populations associated 
with bud outgrowth. We generated two bar-coded sequenc- 
ing libraries using total RNA extracted from sugarcane inac- 
tive buds (0 DAP) and from outgrown axillary buds (2 DAP; 
Fig. 1). After excluding poor-quality reads, those without 
inserts, and those with inserts smaller than 19 nt, we identi- 
fied reads matching known cellular sRNAs of 19-25 nt of 
various classes (Supplementary Table S2 at JXB online). 

The distribution of sequences in different size classes of 
sRNAs was comparable in inactive and developing buds 



(Supplementary Fig. SI A at JXB online), which suggests that 
multiple sRNA pathways are similarly active in cells not under- 
going substantial proliferation (in 0 DAP inactive buds) when 
compared with cells from 2 DAP outgrown buds. As generally 
described, the majority of the sRNA sequences in both librar- 
ies are 24 nt (4 918 279 reads), followed by 21 nt (2 108 665 
reads) and 22 nt (1 478 716 reads) sequences. To understand 
better the possible roles of this complex pool of sRNAs, we 
initially evaluated the distribution of the two major classes of 
axillary bud-associated sRNAs (21 and 24 nt sRNAs) on pub- 
licly available sugarcane BACs that are fully annotated (Jannoo 
et al, 2007; Garsmeur et al, 201 1). We observed a similar pat- 
tern of distribution of these two classes of sRNAs between 
homoeologous haplotypes from Saccharum officinarum and 
Saccharum spontaneum, the ancient polyploid progenitors of 
sugarcane modem commercial hybrids (Dillon et al, 2007). 
Several sequences from the 21 and 24 nt classes of sRNAs 
matched transposable elements (TEs) housed into intergenic 
regions or within protein-coding regions (Supplementary Fig. 
SIB). This indicated that most of these sRNAs are rasiRNAs 
and that these TEs might be transcriptionally active in ancient 
species as well as in modern hybrids. 

Numerous 20-24 nt sRNAs from axillary buds match 
MITEs (Supplementary Fig. SIB), providing evidence for 
active sRNA biogenesis derived from these types of TE in 
the sugarcane genome. MITEs often exist in high copy num- 
ber in genomes and some are transcribed and capable of gen- 
erating sRNAs (Kuang et al, 2009). MITE-derived sRNAs 
may have important roles in sRNA-mediated gene regulation 
and the evolution of several species. Moreover, they might 
represent the evolutionary link between miRNA and siR- 
NAs (Piriyapongsa and Jordan, 2008; Zanca et al, 2010; 
Li et al, 2011). To identify transcribed sugarcane MITEs 
that perfectly matched sRNAs from axillary buds, we first 
searched for MITEs in the EST database (see Materials and 
methods). We identified 23 ESTs similar to known MITEs 
(the majority belonging to the Tourist-like superfamily; data 
not shown), nine of which had perfect matches with axillary 
bud-expressed sRNAs. These MITEs were named as MiScl- 
MiSc9 (Supplementary Table S3 at JXB online), similar to 
the nomenclature provided by Kuang et al, (2009). 
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Fig. 1. Sugarcane axillary bud outgrowth. Sugarcane hybrid SP80-3280 main stem containing inactive or 0 DAP buds (A), and stem 
cuttings containing developing buds at 2 DAP (B) and 5 DAP (C). Tissues from 0 and 2 DAP developing buds were harvested for 
constructing the sRNA libraries, gene-expression analyses, and hormone profiling. Tissues from 5 DAP developing buds were used for 
gene-expression analyses and hormone profiling. Bars, 1 .0cm. (This figure is available in colour at JXB online.) 
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Whilst most transcribed sugarcane MITEs (Tourist-like 
superfamily) had hits with 22-24 nt sRNAs, MiSc8 and MiSc9 
seem to generate mostly 20 and 21 nt sRNAs (Supplementary 
Table S3). MiSc9 displayed a weak similarity to Stowaway- 
like element and formed a harpin-like structure from which 
sRNAs from different sizes were apparently excised at spe- 
cific locations like canonical miRNAs (Fig. 2A). We evalu- 
ated the accumulation pattern of the most represented 21 nt 
sRNA (namely sRNA2908; Supplementary Table S3) that 
initiated with a 5 '-terminal U residue and perfectly matched 
MiSc9 stem-loop via stem-loop RT-PCR (Varkonyi-Gasic 
et al, 2007). sRNA29098 transcripts were readily detect- 
able, indicating that it was an authentic sRNA. Moreover, 
it accumulated in all interrogated monocots but not in the 
eudicots Arabidopsis and tomato, indicating it is poten- 
tially generated only by the genome of monocots (Fig. 2B). 
As sRNA29098 perfectly matched Saccharum spontaneum 
BAC (position around 90 000 bp; Supplementary Fig. SIB), 
we tested whether this sRNA accumulated in Saccharum 
officinarum and Saccharum spontaneum. sRNA29098 tran- 
scripts accumulated in axillary buds from ancient species 
and modern hybrids (Fig. 2B, C). Based upon the fact that 
sRNA29098 and its variants matched a specific location in 
the stem-loop (Fig. 2A), reminiscent of DCL1 cleavage prod- 
ucts, we propose that they originated from MiSc9 via miRNA 
or miRNA-like biogenesis pathways. 

Among the sugarcane ESTs with sRNA hits, it is pos- 
sible that some represent transcripts capable of yielding 21 
nt tasiRNAs (Allen et al., 2005). To investigate the presence 
of this type of siRNA in our axillary bud libraries, we used 
the computational approach described by Chen et al. (2007). 



We identified seven potential tasiRNA-generating ESTs 
(Table 1), two of which represent sugarcane TAS3 precur- 
sors (Shen et al., 2009). We detected a total of nine unique 
bud-expressed tasiRNAs from the two TAS3 precursors, 
one of which (5'-UCUUGACCUUGUAAGACCCAA-3') 
is homologous to Arabidopsis tasiR-ARFs (Allen et al., 
2005) and possibly targets sugarcane AUXIN RESPONSE 
FACTOR (ARF) genes (data not shown). These results sug- 
gested that the TAS3 tasiRNA pathway was operating in sug- 
arcane axillary buds. 

Identification of novel and known miRNAs expressed in 
sugarcane vegetative axillary buds 

Conserved families of miRNAs are found in many plant spe- 
cies and are known to have important functions in plant devel- 
opment (Voinnet, 2009). Conversely, the functions of most 
non-conserved and species-specific miRNAs remain to be 
determined. We identified two new 21 nt miRNA candidates 
(miRcandl and miRcand2) that are generated by distinct 
precursors, both supported by singleton ESTs (Table 2 and 
Supplementary Table S4 at JXB online). Although miRNA* 
accumulation provides strong supporting evidence for the 
cleavage by DCLs during miRNA biogenesis (Voinnet, 2009), 
we could not detect any sequence in our axillary bud libraries 
for the novel miRNA candidates. This observation was not 
surprising as non-conserved miRNAs and their miRNA*s 
are generally expressed at low levels or in specific cell types or 
under specific growth conditions (Rajagopalan et al, 2006). 

BLASTn analysis against all nucleotide sequences in the 
NCBI database revealed that no homologues for miRcand 1 were 
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Fig. 2. MITE-derived sRNAs are expressed in inactive and developing vegetative buds. (A) Predicted stem-loop structure of MiSc9 
shows the SRNA29098 sequence in red. Red lines represent the sRNA reads shown in Supplementary Table S3. (B) Stem-loop pulsed 
RT-PCR to detect sRNA29098 transcripts in tissues of sugarcane [Saccharum spp.), sorghum (Sorghum bicolor), rice (0. sativa), 
brachypodium (B. disthacyon), tomato (Solarium lycopersicum), and Arabidopsis (A. thaliana). The highly conserved miR156 was 
used as an internal control. (C) Stem-loop pulsed RT-PCR of sRNA29098 in inactive (0 DAP) and developing (5 DAP) axillary buds 
of Saccharum officinarum and Saccharum spontaneum. The sugarcane 25S rRNA gene (TC1 48086) was used as a loading control. 
Reactions without RT (-RT) and without cDNA (-) were used as negative controls. 
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Table 1. Putative tasiRNA-generating sugarcane ESTs. 



No. 


EST 


Start 


Strand 


Total no. siRNAs 


No. non-redundant 
tasiRNAs 


P value 


Annotation 


1 


CA151010 


22 


1 


27 


5 


2.67E-04 


Similar to sugarcane TAS3 
precursor EU327139.1 C 


2 


TC1 40837 


449 


1 


21 


5 


3.77E-04 


Similar to sugarcane TAS3 
precursor EU327139.1 


3 


CA1 95273 


305 


-1 


8 


6 


1 .97E-04 


Putative NBS-LRR disease 
resistance protein 






119 


1 


15 


8 


5.52E-05 




4 


CA261958 


407 


1 


80 


11 


1.81E-04 


Putative transposase protein 


5 


DV735856 


216 


-1 


37 


11 


7.75E-04 


Hypothetical protein 






93 


1 


36 


11 


9.45E-04 




6 


DV548691 


69 


1 


91 


15 


6.84E-04 


Saccharum hybrid cultivar 
R570 retrotransposon 


7 


TC1 1 9297 


439 


1 


10 


5 


3.77E-04 


DNA-directed RNA polymerase 
subunit RPABC4-like 



a Total number of 21 -nt tasiRNAs that perfectly match each sugarcane ESTs. 

b Only ESTs perfectly matching at least five non-redundant tasiRNAs (P <0.001) were considered. 

c EST identified by Shen et al. (2009). 



found in other plant species, suggesting that this putative newly 
miRNA is sugarcane specific. In contrast, sugarcane miRcand2 
seemed to have a homologue precursor in the sorghum genome. 
We identified a highly similar sorghum EST (CD227606; 85% 
nucleotide sequence identity) that was analysed further by our 
in-house MIRcheck-based script (Supplementary Fig. S2 at 
JXB online). The sorghum miRcand2 precursor was located 
at chromosome 5 (position 60199113-60199773). To verify 
experimentally whether miRcand2 is generated by other plant 
species, we monitored the accumulation of this 21 nt putative 
novel miRNA in some monocots and eudicots via stem-loop 
RT-PCR (Varkonyi-Gasic et al., 2007). miRcand2 was readily 
detectable by stem-loop RT-PCR, although only in sugarcane 
and sorghum tissues (Fig. 3A). The result corroborated the 
hypothesis that this newly microRNA candidate is restricted to 
the genomes of sugarcane and its closest relatives. Sugarcane 
miRcand2 was expressed in vegetative axillary buds as well as in 
the vegetative apex (apical meristem plus a few leaf primordia) 
of hybrid SP80-3280 at variable levels (Fig. 3B), indicating that 
it may have roles in regulating gene expression during sugarcane 
early developmental stages. 

To gain more insights into the functions of these two novel 
miRNAs in sugarcane, we computationally predicted possi- 
ble targets among the bulk of sugarcane EST sequences (see 
Materials and methods). As a result, predicted target genes 
were identified only for miRcand2, which were represented by 
ESTs encoding pentatricopeptide repeat (PPR)-like proteins 
(Supplementary Table S5 at JXB online). None the less, our 
5' rapid amplification of cDNA ends (5 'RACE) experiments 
for these ESTs were inconclusive (data not shown). Therefore, 
at this point, miRcan2 cannot be confirmed as a functional 
novel miRNA targeting PPR genes in sugarcane. To iden- 
tify known miRNAs in our bud libraries, we employed a 
homology-based approach using BLASTn comparison, miR- 
Base (version 17.0) as a reference set, and sugarcane miRNA 
precursors (Zanca et al, 2010). We found several miRNA 



candidate sequences that exhibited perfect or near-perfect 
matches with at least 26 families of known miRNAs (Table 2 
and Supplementary Table S4). Candidate miRNA reads 
were usually 21 nt and several could be mapped to sugarcane 
miRNA precursors (Supplementary Table S4). In addition to 
the sugarcane miRNA precursors already deposited in miR- 
Base, we identified a novel miR166 precursor that was vali- 
dated by our in-house MIRcheck-based script and matched 
several miR166 candidate reads (Supplementary Fig. S2 
and Table S4). Moreover, 16 known miRNAs identified in 
the developing axillary buds (Table 2) were also expressed 
in leaf tissues of sugarcane plants exposed to drought stress 
(Ferreira et al, 2012). 

Interestingly, we identified highly represented miRNAs in 
our libraries with variants in the 5' or 3' terminus, as well as 
internal mismatches. For instance, the miR159 family had 56 
variants in both libraries (Supplementary Table S4). All poly- 
morphic sequences corresponding to the mature miRl 59 were 
precisely excised at the same specific location of the sugarcane 
MIR159 precursors (Supplementary Fig. S3 at JXB online). 
Although we cannot rule out the presence of sequencing 
errors, the high number of polymorphic sequences is prob- 
ably a reflection of the high level of SsMIR159 gene expres- 
sion in axillary buds combined with the elevated number of 
alleles generally observed for most sugarcane loci (Cordeiro 
et al., 2000). Similar results have been observed for miRNAs 
from allotetraploid cottons (Pang et al., 2009). Whether or 
not these variants have any biological impact on target regu- 
lation in polyploid species remains to be seen. 

Expression profiles of known miRNAs and targets 
during axillary bud outgrowth 

In general, sRNAs associated with TEs and protein-encod- 
ing genes were similarly represented in both inactive and 
developing bud libraries (data not shown). Conversely, most 
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Table 2. Known and novel sugarcane microRNA candidates. 



Number 


miRNA family 3 


No. transcripts (TPM) from 
inactive bud (0 DAP) 


No. transcripts (TPM) from 
active bud (2 DAP) 


Annotation 0 


Known miRNAs 










1 


miR156 


2.5 


5.8 


miRbase, EST 


2 


miR159* 


100244.6 


46877.7 


miRbase, EST 


3 


miR160* 


31.5 


227.8 


miRbase 


4 


miR162* 


88.5 


117.7 


miRbase 


5 


miR164* 


6.1 


34.6 


miRbase 


6 


miR166* 


912.6 


2224.2 


miRbase, EST (CN607727) 


7 


miR167* 


7.8 


11.5 


miRbase, EST 


8 


miR168* 


1157.1 


3937 


miRbase, EST 


9 


miR169* 


6.6 


3.7 


miRbase, EST 


10 


miR171* 


22.4 


450.8 


miRbase 


11 


miR172 


1.2 


3.3 


miRbase 


12 


miR319* 


2584.1 


766.8 


miRbase, EST 


13 


miR393 


7.8 


4.7 


miRbase 


14 


miR394 


8.7 


14.0 


miRbase 


15 


miR395 


3.0 


7.5 


miRbase 


16 


miR396 


91.0 


74.9 


miRbase, EST 


17 


miR397* 


0.0 


6.8 


miRbase 


18 


miR398" 


0.0 


5.4 


miRbase 


19 


miR399 


5.5 


1.6 


miRbase 


20 


miR408* 


0.8 


9.1 


miRbase, EST 


21 


miR444* 


0.5 


6.8 


miRbase, EST 


22 


miR444b.2* 


0.9 


3.5 


miRbase, EST 


23 


miR444d.3* 


2.3 


6.8 


miRbase, EST 


24 


miR529 


18.2 


20.6 


miRbase 


25 


miR827* 


27.6 


150.4 


miRbase, EST 


26 


miR1878* 


39.9 


10.3 


miRbase 


Novel miRNA candidates 










1 


miRcandl* 


30.6 


1.2 


EST (CA257509) d 


2 


miRcand2* 


29.4 


530.6 


EST (CA222602) d 



a microRNAs with an asterisk showed significant differential expression (P <0.01) using the Bonferroni-corrected P value of Fisher's exact test. 
b The normalized expression of transcripts per million (TPM) was obtained by the following formula: (total read count of each miRNA/total 
count of high-quality small RNA reads)x1 000 000. The normalized values were rounded to the nearest tenth. 
c Annotation of the reads based on miRbase versionl 7 and/or matching with EST sequences. 
d Identified in this work. 
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non-conserved miR529 were used as internal controls. Reactions without RT (-RT) and reaction without cDNA (-) were used as negative 
controls. (B) Detection of miRcand2 transcripts through stem-loop pulsed qRT-PCR in 0, 2, and 5 DAP developing axillary buds, as well 
as in the vegetative apex of sugarcane hybrid SP80-3280. All qRT-PCR experiments used inactive bud (0 DAP) as the reference sample 
(set to 1 .0). Error bars indicate standard deviation of three biological replicates. Distinct letters indicate significant differences at P <0.05 
(Student's Mest). R.E., relative expression. 
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miRNAs (20/28) were differentially expressed (P <0.01) dur- 
ing bud outgrowth (Table 2). Post-transcriptional mecha- 
nisms that operate at the biogenesis level could modulate 
the spatiotemporal generation of miRNAs, perhaps through 
differential expression of the necessary biogenesis factors 
(Nogueira et al., 2009). This seems not to be the case for bud- 
associated miRNAs, as we detected similar levels of SsDCLl 
(Saccharum spp DCL1) transcripts in inactive and developing 
buds (Fig. 4A, B). Although we cannot exclude the possibility 
that other miRNA processing components are differentially 
expressed in developing buds, it is likely that the differential 
accumulation of miRNAs in axillary buds is mainly regulated 
at the transcriptional level. 



To complement the sRNA sequencing analysis and assem- 
ble a comprehensive repository of information on miRNAs 
associated with axillary bud outgrowth, we employed stem- 
loop real-time qRT-PCR using total RNA from new biologi- 
cal samples of sugarcane axillary buds. Stem-loop qRT-PCR 
is a standard technique to accurately quantify miRNAs 
(Varkonyi-Gasic et al, 2007). Initially, we compared the 
expression profiles of four sugarcane genes commonly used 
as reference genes in qRT-PCR experiments (Iskandar et ah, 
2004). Based on the qPCR results for three biological rep- 
licates, the 25S rRBA gene (TC 148086) was the most stable 
gene (Supplementary Fig. S4). Therefore, it was selected as 
the reference gene for all qRT-PCR experiments in this study. 
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Fig. 4. DCL1 and known miRNAs are expressed during axillary bud outgrowth. (A) Detection of SsDCLl transcripts through qRT-PCR 
in 0, 2, and 5 DAP developing axillary buds of sugarcane hybrid SP80-3280. (B) Inferred phylogenetic relationships among SsDCLl 
orthologues in rice (Os) and Arabidopsis (At). The p-distances were calculated from amino acid alignments of conserved blocks, and 
tree topology was inferred with the neighbour-joining method. Only bootstrap values higher than 50% are shown for 1000 replicates. 
SsDCLl is underlined. (C) Accumulation of selected mature miRNAs was quantified through stem-loop pulsed qRT-PCR in 0, 2, and 5 
DAP axillary buds of sugarcane hybrid SP80-3280. All qRT-PCR experiments used inactive bud (0 DAP) as the reference sample (set 
to 1 .0). Error bars indicate standard deviation of three biological replicates. Distinct letters indicate significant differences at P <0.05 
(Student's f-test). R.E., relative expression. 
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To verify the expression profiles of miRNAs in inactive and 
developing buds (2 and 5 DAP buds), we selected six broadly 
conserved and two non-conserved miRNAs for further 
analysis. All were readily detected by stem-loop qRT-PCR 
(Fig. 4C), which indicated that these miRNAs are authentic 
sRNAs. Whilst miR159, miR319, and miR444 transcripts 
accumulated at higher levels in inactive buds, miR156 and 
miR171 were expressed at higher levels in outgrown buds. 
miR529 and miR160 transcripts accumulated at low levels 
specifically in 5 DAP developing buds (Fig. 4C). Some dis- 
crepancies between the frequency of miRNAs represented 
in the libraries and the stem-loop qRT-PCR data (Table 2 
and Fig. 4C) might have resulted from insufficient depth of 
sequencing coverage or cloning bias. 



In this study, we computationally predicted more than 30 tar- 
gets for sugarcane miRNAs. All sugarcane cDNAs predicted to 
be targets of known miRNAs were homologues of miRNA tar- 
get genes in Arabidopsis and other plant species and may be asso- 
ciated with various developmental processes (Supplementary 
Table S5). We monitored the expression of five homologous 
miRNA targets by qRT-PCR: Saccharum spp. SQUAMOSA 
PROMOTER BINDING-LIKE PROTEIN1 (SsSPLI), 
SsGAMYB, SsTCPl, AUXIN RESPONSE FACTOR10I16 
(SsARFI0/16), and SCARECROW-LIKE1 (SsSCLl) 
(Fig. 5A). Interestingly, all the genes were upregulated at varia- 
ble levels in developing buds when compared with inactive buds 
(Fig. 5B). Among them, SsGAMYB displayed an inverse rela- 
tionship with miR159 expression at all developmental stages 
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Fig. 5. Expression profiles of selected miRNA targets during axillary bud outgrowth. (A) miRNA-complementary sites in the target 
mRNA and the miRNA. Watson-Crick pairing (vertical dashes) is indicated. (B) Detection of SsSPLI, SsGAMYB, SsTCPl, SsARF10/16, 
and SsSCU transcripts was performed by qRT-PCR in 0, 2, and 5 DAP axillary buds of sugarcane hybrid SP80-3280. All qRT-PCR 
experiments used inactive bud (0 DAP) as the reference sample (set to 1 .0). Error bars indicate standard deviation of three biological 
replicates. Distinct letters indicate significant differences at P <0.05 (Student's f-test). R.E., relative expression. 
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(Fig. 4C and Fig. 5B). In Arabidopsis, the levels of miR159 and 
GAMYB transcripts are oppositely regulated by the phytor- 
mone ABA (Reyes and Chua, 2007). To examine whether there 
was a correlation between ABA levels and miR159/ SsGAMYB 
expression profiles, we quantified ABA and catabolites in sug- 
arcane axillary buds. ABA levels were reduced from the inac- 
tive buds to 2 and 5 DAP buds (Supplementary Fig. S5). The 
sharp increase in the levels of dihydrophaseic acid suggested 
that bioactive ABA was initially biosynthesized in inactive buds 
and further catabolized in developing buds mainly through the 
8'-hydroxylation pathway (which results in phaseic acid that is 
further reduced to dihydrophaseic acid). These data suggested 
that ABA content within axillary buds correlates positively with 
miR159 accumulation, whilst the opposite was observed for 
SsGAMYB transcript levels. 

Expression patterns of miRNAs in sugarcane 
axillary buds 

To gain further insights into the possible roles of sRNAs 
during sugarcane axillary bud development, we investigated 



the spatial and temporal patterns of miRNA accumulation 
by in situ hybridization. We selected two microRNAs for our 
analysis: a highly represented miRNA (miR159) and a low- 
level detected microRNA (miR156) based on our axillary bud 
sRNA libraries (Table 2). Consistently with its abundance in 
the sRNA libraries, miR159 display high expression levels in 
axillary buds (Fig. 6A, B). miR156 was broadly expressed in 
axillary meristems and leaf primordia at low levels (Fig. 6C), 
whilst miR159 showed high expression levels in vascular bun- 
dles of young leaf primordia (Fig. 6B). The expression pat- 
tern of miR159 was similar in inactive buds (Fig. 6A) and 5 
DAP developing buds (Fig. 6B), suggesting that its transcript 
levels rather than spatiotemporal patterning were develop- 
mentally regulated during axillary bud outgrowth. 

Discussion 

Axillary bud outgrowth is crucial to control plant shoot 
architecture, which is important for biomass production of 
biofuel crops such as sugarcane. Although it is well known 
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Fig. 6. Spatiotemporal expression patterns of miRNAs in axillary buds of sugarcane. (A, B) A probe of a 3'-labelled LNA-modified 
oligonucleotide detecting miR159 was hybridized with longitudinal sections of 0 DAP (A) and 5 DAP (B) vegetative buds. (C) A similar 
probe detecting miR156 was hybridized with longitudinal sections of 0 DAP buds. (D) A scramble-miRNA 3'-labelled LNA probe was 
used as a negative control. Purple and brown staining show probe localization. Arrows represents meristem positions and arrowheads 
indicates young leaf primordia. Bars, 10 urn. 
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that bud outgrowth is controlled mainly by conserved net- 
works that comprise a complex set of transcription factors 
(Leyser, 2009), the presence of additional layers of gene reg- 
ulation embedded in these networks remains poorly under- 
stood. Therefore, the identification of sRNAs and elucidation 
of their functions in vegetative bud development will help us 
to understand better the genetic networks underlying plant 
shoot architecture. 

To identify novel genetic factors that potentially play roles 
in vegetative bud outgrowth, we employed a deep-sequencing 
approach coupled with conventional RNA detection meth- 
ods to profile sRNAs in inactive and developing axillary buds 
of the commercial sugarcane hybrid SP80-3280. We have pro- 
vided the first draft of the sRNA transcriptome in sugarcane 
vegetative axillary buds, which reinforces the idea that the 
sRNA component is complex and diverse within the plant 
kingdom. Although high-throughput sRNA profiling is gen- 
erally straightforward for model plants (Zhu et al, 2008), it 
is rather challenging for plant species such as the polyploid 
sugarcane for which the genome is highly complex and largely 
unknown. However, even using limited publicly available 
genomic and EST sources, we were able to identify known 
and novel miRNAs and show that several siRNAs expressed 
in buds match distinct TEs and that some are phased tasiR- 
NAs originating from specific transcripts. 

Diverse populations ofsRNAs are expressed in inactive 
and developing axillary buds 

In general, we found similar sRNA abundance when com- 
paring inactive and growing axillary buds, suggesting that 
inactive buds are also actively generating sRNAs to control 
gene expression. This data is consistent with the idea that 
quiescent cells in inactive/dormant buds are metabolically 
and molecularly active (Stafstrom et al, 1998). For instance, 
several bud-expressed rasiRNAs matched intergenic and 
intron-embedded TEs, indicating that these TEs are either 
active or under RNA-dependent epigenetic regulation. Such 
epigenetic regulation may be crucial for the developmental 
programmes associated with outgrowth of vegetative axillary 
buds, similarly to chromatin modifications that take place 
during bud dormancy release in some woody perennial plants 
(Leida et a!., 2012). 

MITEs, a type of non-autonomous TE, may assert gene 
regulation via sRNA pathways (Kuang et al, 2009). Our 
identification of MiSc sRNAs indicated that MITEs could 
serve as both origins of sRNAs and targets for sRNA-medi- 
ated gene silencing during axillary bud outgrowth of mod- 
ern sugarcane hybrids and ancient progenitor species. For 
instance, sugarcane genes bearing MITE insertions may 
be co-opted for using sRNA-based regulation, similar to 
regulation imposed by SINE-like direct repeats in the FWA 
promoter (Chan et al, 2004). Additionally, some MiSc rasiR- 
NAs expressed in axillary buds are probably generated and 
may act similarly to miRNAs (Fig. 2), as has been shown 
recently for two rice MITE-derived sRNAs that positively 
regulate ABA signalling by post-transcriptionally regulating 
MAIF1 expression (Yan et al, 2011). All these mechanisms 



are expected to contribute to the complexity of epigenetic 
regulation during axillary bud outgrowth. 

Four families comprising eight tasiRNA loci have been 
described in Arabidopsis, whereas hundreds of loci of 
unknown function generate phased sRNAs in grasses 
(Johnson et al, 2009). In sugarcane axillary buds, we identi- 
fied tasiRNAs that seemed to be generated by a small number 
of transcripts (Table 1), including TAS3 and nucleotide-bind- 
ing site-leucine-rich repeat (NBS-LRR) disease resistance 
protein-coding transcripts. Secondary tasiRNAs generated 
from these types of transcript are known to regulate gene 
expression in trans in several plant species (Allen et al, 2005; 
Zhai et al., 2011). Interestingly, we detected sequences similar 
to 21 nt miR390 (sRNA that triggers tasi-ARF production; 
Allen et al, 2005) and 22 nt miR2118 (sRNA that triggers 
the production of NBS-LRR tasiRNAs; Zhai et al, 2011). 
However, as there were fewer than 15 reads representing these 
miRNAs in our sRNA libraries, they were not included in 
further analyses (see Materials and methods). The tasi-ARF 
pathway has been shown to play a role in leaf development 
(Nogueira et al, 2007), whilst NBS-LRR tasiRNAs are 
involved in pathogen defence pathways in leguminous plants 
(Zhai et al, 2011). Our data suggest that sugarcane TAS3 
and NBS-LRR gene-derived tasiRNAs are also part of the 
intricate regulatory networks operating during vegetative bud 
development. It will be interesting to test whether siRNAs 
from the additional transcripts (Table 1) act as trans-acting or 
cw-acting siRNAs in sugarcane. 

Dynamic expression of miRNAs and targets during 
axillary vegetative bud outgrowth 

Comparison of expression profiles between miRNAs and the 
corresponding target genes can reveal positive or negative 
correlations, which expand our understanding of how miR- 
NAs are involved in sugarcane vegetative bud development. 
Specific microRNAs displayed dynamic expression profiles in 
inactive and developing vegetative buds (Fig. 4C), suggest- 
ing that they cooperate with their targets to achieve global 
gene-expression networks that function within axillary buds. 
Interestingly, the expression profiles of some miRNAs were 
positively correlated with their targets. For instance, levels of 
transcripts of miR171 and one of its potential targets in sug- 
arcane, SsSCLl, increased as the bud grew out (Figs 4C and 
5B). This may due to the consequence of feedback regulation 
(Baulcombe, 2004) or to the fact that this miRNA may func- 
tion by translational repression (Lanet et al, 2009). 

SsSPLl is a member of the SQUAMOSA PROMOTER 
BINDING-LIKE PROTEIN (SPL) family of transcription 
factors, of which some members are associated with shoot 
development and plant architecture (Schwarz et al, 2008; 
Jiao et al, 2010). In monocots, some SPL mRNAs are tar- 
geted by both miR156 and miR529 (Fig. 5A; Chuck et al, 
2010). Based on our qRT-PCR data, levels of miR529 and 
SsSPLl transcripts were inversely correlated specifically in 5 
DAP developing buds, consistent with SsSPLl being nega- 
tively regulated by this miRNA. Interestingly, the expres- 
sion profile of miR156 was positively correlated with that 
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of SsSPLl at all developmental stages (Figs 4C and 5B). 
Moreover, miRl 56 was less represented in both libraries when 
compared with miR529 (Table 2). It is possible that SsSPLl 
expression in axillary buds is independent of both miRNAs 
in the early stages of development. However, at the late stages 
of bud outgrowth (5 DAP), SsSPLl transcript levels may be 
regulated predominantly by miR529. 

miRNA159 was the most enriched miRNA in our axillary 
bud libraries, followed by miR319 (Table 2). The specificity of 
these closely related miRNAs is achieved in part by their dis- 
tinct expression, with miR319 expressed at much lower levels 
than miRl 59 in Arabidopsis tissues (Palatnik et al, 2007). 
Consistently, miR319 was underrepresented in our sRNAs 
libraries when compared with miRl 59 (Table 2). Interestingly, 
miRNA159 and miR319 accumulated at high levels in inac- 
tive buds (Table 2), whilst their targets (SsGAMYB and 
SsTCPl, respectively) were upregulated in developing buds 
(Figs 4C and 5). 

In situ localization of miRl 59 in axillary buds revealed 
that it accumulated mainly in young leaf primordia (Fig. 6A, 
B) and could potentially regulate its targets in these tissues. 
miR159-regulated GAMYB-Yks genes encode R2R3 MYB 
domain transcription factors that are implicated in ABA 
signalling in germinating seeds (Reyes and Chua, 2007). 
Interestingly, the ABA content of vegetative axillary buds is 
closely correlated with bud dormancy/inactivity (Shimizu- 
Sato and Mori, 2001). Although speculative, we propose a 
scenario in which miRl 59 levels in sugarcane axillary buds 
are positively regulated by ABA, as has been shown in ger- 
minating seeds of Arabidopsis (Reyes and Chua, 2007). The 
differential accumulation of ABA in inactive and outgrown 
buds (Supplementary Fig. S5) may converge into a fine-tuned 
regulation of miRl 59 expression, which, in turn, would act 
as a molecular switch during the inactive stage of bud devel- 
opment, permitting higher expression of SsGAMYB only in 
growing axillary buds. Further functional studies are needed 
to clarify the roles of miR159-regulated pathways associated 
with axillary bud development. 

In summary, we have presented a comprehensive survey 
of sRNAs in vegetative axillary buds of the tropical biofuel 
crop sugarcane. Collectively, the dynamic changes in miRNA 
accumulation and their targets probably promote signalling 
and metabolic pathways that are essential for bud outgrowth. 
Future research on the miRl 59 family and target regulation 
may increase our understanding of the biological role of this 
miRNA and its contribution to the early emergence of vegeta- 
tive buds. Our work has opened a new avenue for functional 
studies on sRNA-mediated gene regulation in vegetative 
axillary bud outgrowth and consequently plant vegetative 
architecture. 
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to MITEs and their associated sRNA sequences. 

Supplementary Table S4. Total reads of known and novel 
sugarcane miRNA candidates identified in the axillary bud 
libraries. 

Supplementary Table S5. Predicted target sequences of 
known and novel sugarcane miRNAs identified in the axillary 
bud libraries. 
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